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Abstract 

We present a Ktt mass spectrum analysis of the four-body semileptonic charm 
decay K~TT^pAv in the range of 0.65 GeV/c^ < mKvr < 1-5 GeV/c^. We 

observe a non-resonant contribution of 5.30 ± 0.74lQgg% with respect to the to¬ 
tal D~^ —> K~tt^ decay. For the A*(892)^ resonance, we obtain a mass of 
895.41 ± 0.32 1q43 MeV/c^, a width of 47.79 ± 0.86lj;'Qg MeV/c^, and a Blatt- 
Weisskopf damping factor parameter of 3.96 ± 0.54Iggg GeV“^. We also report 

90% CL upper limits of 4% and 0.64% for the branching ratios ^ 
and ^^r(D+^K-^n+\+u)'^ ’ respectively. 
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Weak semileptonic decays of charm mesons continue to attract interest due 
to the relative simplicity of their theoretical description: the matrix element 
of these decays can be factorized as the product of the leptonic and hadronic 
currents. This makes the ^ K~t [^decay a natural place to study the 
Ktt system in the absence of interactions with other hadrons. Due to Watson’s 
Theorem [1,2], the observed Kn phase shifts in should be 

the same as those measured in Kti elastic scattering. 

It is known that the Ktt hnal state of decay is strongly dom¬ 

inated by the i^*(892)° vector resonance [3,4]. The large and clean sample of 
D~^ —events collected by the Fermilab FOCUS experiment pro¬ 
vides an excellent opportunity to measure the iF*(892)° mass and width, as 
well as the effective Blatt-Weisskopf damping factor parameter discussed in 
Ref. [5]. We also search for structures other than the W*(892)° resonance in 
the mass range of 0.65 GeV/c^ < niKit < 1-5 GeV/c^. 

The hrst suggestion that the ^ K~ti^ fi'^v decay proceeds via states other 

than the iF*(892)° resonance comes from the Fermilab E687 experiment [4]. 
The presence of an additional structure was conhrmed by FOCUS in the anal¬ 
ysis of the angular decay distributions, in which the iF*(892)° form factor 
was measured [6,7]. Specihcally, signihcant discrepancies were found between 
the data and the predicted iF*(892)V"'"^ angular decay distributions. 

A nearly constant amplitude and phase contribution to the helicity zero am¬ 
plitude of the virtual was required to adequately £t the observed de¬ 
cay angular distributions. The s-wave amplitude, 006**^°, was measured in the 
vicinity of the iF*(892)° pole with parameters Oq = 0.330 ± 0.022 ± 0.015 
and (5o = 0.68 ± 0.07 ± 0.05. This new component accounts for 5% of the 
^ branching fraction. 

Motivated by this earlier FOCUS result, we search for other contributions 
in the Kti spectra. Specifically, we look for a possible contribution from the 
77*(1680)°, i7g(1430)°, and k. We also present a more complete description 
of the non-resonant contribution. The existence of the k, reported in [8], re¬ 
mains controversial due to difficulties in the theoretical treatment of broad 
scalar states and the absence of a clear observation of this state in scattering 
experiments. Many models predicting the decay width of semileptonic decays, 
such as ISGW2 [9] and QCD Sum-Rules [10], indicate the tendency for these 
decays to proceed via low mass structures. In [11] it is suggested that if the 
K, has a substantial qq component in its wave function, it could account for 
more than 10% of the ^ K~n ~^decay rate. 

The data were collected in the Wideband photoproduction experiment FO¬ 
CUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, a forward 
multi-particle spectrometer is used to measure the interactions of high energy 
photons on a segmented BeO target. The FOCUS detector is a large aperture, 
fixed-target spectrometer with excellent vertexing and particle identification. 
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The FOCUS beamline [12] and detector [6,13,14,15] have been described else¬ 
where. 

To isolate D~^ ^ K~Tr~^ events, we require that the muon, pion, and kaon 
candidate tracks have a 5% or greater conhdence level to originate from a 
common secondary vertex. Background is reduced by requiring the secondary 
decay vertex be separated from the production (primary) vertex by greater 
than 10 (Te, where ae is the uncertainty on the separation between the primary 
and secondary vertices. Possible backgrounds from higher multiplicity charm 
decays are suppressed by requiring the K~ti ^vertex be isolated from other 
tracks in the event (excluding tracks from the primary vertex). Specihcally, we 
require that the maximum conhdence level for another track to form a vertex 
with the secondary vertex candidate be less than 1%. To suppress background 
from secondary interactions the decay vertex candidate must he outside any 
target foil or detector material. 

The muon, pion, and kaon candidates are selected in the following way. The 
muon track must have hits in at least 5 of the 6 segmented scintillator layers 
which comprise the inner muon detector and a muon conhdence level exceeding 
5% (based on the ht to the hits). The pion and kaon tracks must have a muon 
conhdence level less than 0.1%. The kaon is required to have a Cerenkov light 
pattern more consistent with that of a kaon than that of a pion by 2 units of 
log-likelihood, while the pion track is required to have a light pattern favoring 
the pion hypothesis over that of the kaon by 2 units [14]. In addition, the 
pions and kaons are required to have momenta greater than 5 GeV/c, while 
the muon momentum must exceed 10 GeV/c. 

To suppress background from K~7r~^7r~^, we require that the invari¬ 

ant mass of the three tracks, where the muon candidate is given the pion 
mass, is less than 1.8 GeV/c^. To suppress background from 
{K ~and D*~^ {K*~, we require m{K~— 

> 0.2 GeV/c^. A total of 18245 candidates re¬ 

main after the selection criteria. 

The charm background, charm decays that are not ^ K~n ~^is esti¬ 
mated from more than one billion charm Monte Garlo events that pass through 
the entire data analysis chain. Our Monte Garlo is based on Pythia [16] and 
incorporates all known charm decays. The charm Monte Garlo sample was 
scaled to the data sample size using the htted yield of the D~^ ^ sig¬ 

nal. To estimate the background contribution coming from non-charm events, 
we dehne a wrong sign sample (WS) formed by K^t [~tracks in the sec¬ 
ondary vertex. We assume that non-charm events populate the wrong sign and 
right sign (RS) samples equally. The non-charm background distribution is ob¬ 
tained by subtracting the WS charm background (obtained from the Monte 
Garlo sample) from the WS data sample. We estimate the charm and non¬ 
charm background contributions to be, respectively, 17.8% and 3.2% of the 
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total number of events over our signal region. 


Four-body decays of spinless particles are described by five kinematic variables. 
The variables chosen in this analysis are the K~tt^ invariant mass {mKn), the 
square of the mass (g^), and three decay angles: the angle between the tt’*' 
and the D~^ direction in the K~7r~^ rest frame {0^), which dehnes one decay 
plane, the angle between the u and the direction in the rest frame 
(6*^), which dehnes the second decay plane, and the acoplanarity angle (y) 
between these two decay planes. 

The differential decay rate can be represented by a coherent sum of resonant 
and non-resonant contributions to the angular momentum eigenstates of the 
K~ti^ system. 


dr 

dfrixTr 


-Mj A. 


MR) 


J R 


(p dVt 


( 1 ) 


where dVL = dq^ dcos 6^ dcos 9^ dx, Af j is the weak matrix element for a transi¬ 
tion with angular momentum J, Aj^i^r) represents the form of the hadronic hnal 
state amplitude contribution of resonance R (or non-resonant) with strength 
Rj,(r)) and (p is the phase space density. 

The possible resonant states that couple to K~7r~^ are the scalars k and 
the vectors K*{892)^ and iF*(1680)°, and the tensor 
The non-resonant contribution is assumed to be scalar. ^ Small amplitude con¬ 
tributions are most likely to be observed through the interference with large 
amplitude components. Due to the orthogonality of states with different an¬ 
gular momentum, only amplitudes with the same spin will produce signihcant 
interference contributions to the mK-w mass spectrum, given our reasonably 
uniform angular acceptance. Therefore, the small vector iF*(1680)° and scalar 
7 ^ 0 ( 1430 )° contributions might produce an observable effect on the ttikit spec¬ 
trum through their interference with the iF*(892)° and a low mass s-wave 
amplitude, respectively. By contrast, the inclusion of a small 7^2(1430)° res¬ 
onance contribution is unlikely to be observed, since it is orthogonal to the 
(dominant) iF*(892)° and low mass s-wave amplitudes. For this reason we 
do not include the 7^2(1430)° resonance in our hts to the ttikti spectrum in 

The parametrization of resonant states with angular momentum J is given by 


^ Although the simplest way to obtain the forward-backward asymmetry described 
in Ref. [6] is to assume an s-wave amplitude interfering with the AT*(892)^ (as was 
done in Ref [6]), small spin 2 components cannot be excluded. 


5 



the product of a Breit-Wigner ^ and the normalized R ^ K tt’*' coupling, jFj 


•Aj^r — 


mo Be 


-ml + imo T{mK-K) 




( 2 ) 


where = Bq ^ P* is the magnitude of the kaon momentum in 

the resonance rest frame, pX = p*(mo), To = 1, and T\ = ^ ■ B is the 

_Po 

Blatt-Weisskopf damping factor given by -B = 1/+ Tq [5]. The damping 
factor adds an additional ht parameter, tq, in our hts to the iB*(892)'^ line 
shape. The line shape of the k resonance is expected to deviate signihcantly 
from a pure Breit-Wigner, due to its large width and the close vicinity of the 
Kn threshold. In this analysis we use the k, line-shape adopted by E791 [8]. 


We use an empirical parametrization from K~n~^ elastic scattering experi¬ 
ments for the non-resonant amplitude. A partial wave analysis performed by 
LASS observed that the s-wave amplitude can be represented as the sum of a 
Aq( 1430)'^ resonance coupled to K~7r~^ and Kp' and a smooth shape, consis¬ 
tent with the non-resonant hypothesis [17]. ^ LASS fitted the non-resonant 
component to an effective range model of the form 

1 bv* 

cot (5 lass = —^ 

ap* I 

where a = 4.03 ± 1.72 ± 0.06 GeV-^ and b = 1.29 ± 0.63 ± 0.67 GeV-B 
Removing the two-body phase space factor, given by ——, from LASS non- 
resonant amplitude, which is already included in Eq. 1, we obtain the following 
parametrization for the non-resonant hadronic hnal state interaction: 

BInr = sin((5LAss)e*'^^^®L (4) 

p* 


The weak matrix element for the vector process, Afi, and for the scalar pro¬ 
cess, Aio, are written as a function of helicity amplitudes. Hi, derived in [18]. 
Neglecting the mass of the charged lepton the matrix elements are 


Ml 



(1 -f cos6 *£) sinH+{q^,mKTr) 


— {1 — cosOi) sin9,,e — 2sm9iCos9^Ho{q‘^,mKTr) 


( 5 ) 


^ The Breit-Wigner form used by FOCUS differs by a factor of —1 from the 
LASS [17] form. 

^ Charm decays are traditionally fit to a model where the strengths of both resonant 
and non-resonant contributions are fit parameters. Hence we will independently 
adjust the non-resonant and Ag(1430)^ resonant contributions found by LASS to 
best fit our data. 


6 



and 


( 6 ) 


Mo = -2 



sin 6*^ H^{q‘^,mK 7 T)- 


The three form factors for the vector states and the one for the scalar states 
are written assuming the single pole dominance ansatz given by: 


funsutzi^Q ) 


/(o) 


i-9VyL=' 


( 7 ) 


The vector states use the nominal spectroscopic pole masses, Ma = 2.5 GeV/c^ 
and My = 2.1 GeV/c^, and the recent form factor measurements in Ref. [7]. 
The scalar states use My = 2.1 GeV/c^ and the respective zero recoil form 
factor is arbitrarily set to one since its value can always be absorbed in the 
amplitude parameter ao,(ij). 

Next we discuss the angular distribution described by Eq. 1. The Kn spec¬ 
trum described by this equation includes the dominant contribution from the 
i7*(892)° resonance, possible high mass contributions from the R’o(1430)'^ and 
77* (1680)° resonances, and low mass scalar components comprised of a non¬ 
resonant and a possible k contributions, both populating the region where 
relevant discrepancies were found. As discussed in [6,7], the ttixyt distribution 
weighted by cos 9^ provides information on the phase of the additional struc¬ 
ture relative to that of the 77* (892)°. It can be used to discriminate different 
combinations of low mass states, given the large difference between their phase 
shifts. ^ Figure 1 compares the distribution obtained in the data with the pre¬ 
dictions from the non-resonant and k models in the absence of additional phase 
shifts. 

Since a simulation using the LASS parametrization of the non-resonant contri¬ 
bution is sufficient to reproduce the data, we exclude a possible k contribution 
from further consideration. 

Having excluded the k, the most general differential decay rate for > 

in uik-k is given by Eq. 8 where e represents the detector acceptance 
and efficiency: 


dr 

dTflXn 


= I elAfiV + Mo^^l^ 0 dn 


( 8 ) 


^ The expected ruxTr distribution weighted by cos 9 ^ for a pure 71+ —>77*(892)°;U+jv 
decay would be nearly zero. 
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Fig. 1. The background subtracted distribution of rriK-K weighted by cos0v The 
data (squares) show good agreement with the LASS non-resonant parameterization 
(solid histogram) but not with a k model (dashed histogram) 

with vector and scalar amplitudes given by 

V = aii-*(892)0‘^iC*(892)0 + Ojc* (1680)0(1680)0 ( 9 ) 

S = OnR-^NR + Or* (1430)0(1430)0- (10) 


The amplitude coefficients ur*( 392)0 ^ nR'*(i 680 ) 0 ) ®nR) and ur"*( 1430)0 are real, as 
required by Watson’s Theorem [2], 

Equation 8 can be conveniently factorized as: 


dr 

d-TflXir 


\V\^Fn + \S\^Foo + 2^{V* S)Foi 


( 11 ) 


where Fjf ^ I </)dQ, are real functions® that depend only on rriKTr- 

The Fjj' functions are computed from the niKTr spectrum obtained from a 
complete simulation of events, generated according to phase 

space and weighted by M}Mj' and thus represent the intensity modified by 
acceptance and efficiency. The three Fjj' functions are shown in Figure 2. The 
I Vp, |iSp, and 3?(V*i5) functions depend on rriK-K as well as on all £t parameters. 
The cross-term, 2 3fJ(V* iS)Toi, represents the interference between the vector 
and scalar contributions. 


® Because of Eq. 5, all imaginary pieces of will appear as sinusoidal func¬ 

tions of y. Hence any imaginary terms vanish when averaged over y given our nearly 
uniform acceptance in this variable. 
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^Kn GeV/c^ 


Fig. 2. Distributions of the relevant Fjji functions which represent the intensity 
modified by acceptance and efficiency and include the effects of phase space and the 
weak matrix elements. The Fn, Fqq, and Fqi distributions are shown as triangles, 
squares, and circles, respectively. 


The contribution from each decay mode, as well as the i^*(892)'^ parameters, 
are obtained from an unbinned maximum-likelihood £t. We dehne the proba¬ 
bility density function as the sum of the probability density for the signal, £ 5 , 
and for the background, £b- The signal density is described by Eq. 11. The 
background density is given by the sum of charm and non-charm contributions. 
The relative contribution of the two background sources as well as the rela¬ 
tive fraction of the background with respect to the selected L>~'~ K~ti^ 
sample, fs, are fixed at the estimated values, described previously. We fit the 
data by minimizing the quantity u, 


0 ; = — 21 n ^ [(1 —/b)£s +/b£b] ( 12 ) 

events 


The fit parameters are the magnitudes of each amplitude in the signal prob¬ 
ability density function (a*), the nominal mass and width of the it'*(892)°, 
and the parameter tq of the Blatt-Weisskopf damping factor. The it'*(892)° is 
taken as the reference amplitude { 892)0 = 1). The parameters of all other 
resonances are fixed to the PDG values [19]. Decay fractions are obtained in¬ 
tegrating each individual amplitude over the phase space and dividing by the 
integral over the phase space of the overall amplitude. 

To account for momentum resolution effects on the iP*(892)° parameters, we 
refit the data fixing all parameters except the iP*(892)° width and use the 
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probability density function, , given by Eq. 13: 


^^(rriK^) ^ f - niK^, a) 


( 13 ) 


The new probability density function, C^{mK-w), represents the convolution of 
the data £t function with a Gaussian distribution, G, with a = 5.88 MeV/c^, 
value obtained from Monte Carlo simulation. The smearing due to momentum 
resolution increases the i^*(892)° width by approximately 2 MeV/c^. 

Using the procedure described above, we fit the data assuming only a 
7^*(892)°/i^z/ process. The confidence level of this fit is 0.21%, indicating the 
need for additional contributions in the decay. 


The inclusion of a non-resonant scalar component, referred to as the NR 
model, significantly improves the confidence level of the fit to 66%. We find 
^if*( 892 )o = 895.41 ± 0.32 MeV/c^, £^*( 892 )° = 47.79 ± 0.86 MeV/G, tq = 
3.96 ± 0.54 GeV“^, onr = 0.327 ± 0.024, which correspond to a scalar frac¬ 
tion of 5.30 ± 0.74 %. Figure 3 illustrates the contribution of both the ^ 
77*(892)°/i+z/ and non-resonant s-wave process to the observed ttirtt spectrum. 


We also consider possible D~^ 77*(1680)°/i’'"^ and D~^ 77Q(1430)°/i’'"^ 

contributions to our model. Since the data is already well described by a 
model having only the 77*(892)° and non-resonant components, we do not 


expect large contributions from these modes. Including both decays we find 
^K*{ 892 )o = 895.0 ± 1.1 MeV/c^, r 7 f*( 892 )o = 47.63 ± 0.91 MeV/c^, tq = 
5.7±4.8 GeV“^, onr = 0.287±0.073, a;^*(ie8o)o ~ —0.16±0.36, and = 

—0.048 ± 0.19. The 77*(1680)° and 77g(1430)° amplitudes are consistent with 


zero and we find + 

r (D+—~ TT+/^+h*) 


< 4.0% and < 0-64% at 


90% GL. The upper limits are calculated using the method described in [20] 
and assume BR(77*(1680)° —= 0.258 and BR(77 q( 1430)°—= 
0.62 [19]. When the 77*(1680)° is included, we observe a strong correlation 


between rg and ^^‘(iggg^o, inflating the errors on both quantities. To study 
the statistical significance of these new amplitudes, we use a hypothesis test 


based on the maximum-likelihood ratio method [21]. This method compares 
two hypotheses and points out unnecessary degrees of freedom. As a result, 
we obtain a confidence level of 80% in favor of the simple NR model. 


We consider several sources of systematic errors. These include variations of 
the fit conditions, split sample errors, and the uncertainty on the absolute mass 
scale of the experiment, relevant for the 77* (892)° mass measurement. Twenty- 
seven variations of the fit procedure are considered. Starting from the final 
sample we adopt more stringent selection criteria changing the significance of 
the separation between secondary and primary vertices, the secondary vertex 
isolation requirement, and the cut on the muon confidence level. In addition, 
we vary the relative fractions of the different background components. We 
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Table 1 

Estimated systematic uncertainty obtained for each component. 


w-x*(892)0 (MeV/c^) 
Tif*( 892)0 (MeV/c^) 
ro (GeV-i) 
Scalar fraction (%) 


Cut and Model Vars. 

Split Sample 

Mass Scale 

+ 0.13 
-0.29 

+ 0.81 
-0.18 

+ 1.16 
-0.67 

+ 0.95 
-0.92 

±0.11 

±1.05 

±0.61 

±0.28 

±0.30 


vary by ila the values of the parameters from the LASS effective range 
parametrization (Eq. 3). We also include as a systematic error the difference 
between the FOCUS results, obtained from the NR model, and the results 
obtained from the model with high mass structures. Errors from this source are 
asymmetrical: we take the difference between the central and highest/lowest 
values of each £t parameter and scaled by 0.68 to obtain the contribution to 
the systematic error. 

The split sample component takes into account systematic effects introduced 
by residual differences between data and Monte Carlo. This component is 
determined by splitting the data into five pairs of independent subsamples, 
according to the charge, data taking conditions, primary vertex multiplic¬ 
ity, muon momentum, and the momentum of the Kti system. The treatment 
used for the split sample is known as Unconstrained Averaging, described 
in [19, page 14]. 

The total systematic error is given by the sum in quadrature of the uncer¬ 
tainties from the independent sources. Table 1 presents the results of the 
systematic uncertainty evaluation for the measurements. 

Table 2 summarizes the results obtained from the fits using the two mod¬ 
els. The values of the fit parameters are compared to the world average val¬ 
ues [4,7,17,19]. Our measurements of the iL*(892)° mass and width are both 
more than 1 a below the PDG average values. Figure 4 shows a comparison 
between our standard NR model with free K*{892)^ parameters and a NR 
model with the mass and width of 77* (892)° resonance fixed to the world av¬ 
erage values [19]. With the inclusion of a non-resonant contribution, the value 
we obtain for the Blatt-Weisskopf parameter is consistent with LASS [17]. 
The fraction of the scalar component is compatible with the value obtained 
previously in the analysis of the cos 6^ asymmetry [7]. 

In conclusion we have measured the 77*(892)° parameters using a large sample 
of Tl'*'—77*(892)°/i+z/ signal events over a wide mass range. The absence of 
high mass resonances as well as the small background contribution provides 
a unique environment to study the 77*(892)° mass and width. The 77*(892)° 
mass and width measurements are stable with respect to model variation. Our 
measurements of the mass and width are more than 1 a below the present world 
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Fig. 3. Fit to the mK-w data using the NR model. The error bars, the solid lines, the 
dashed lines, and the dotted lines correspond to the data, the model, the background 
contribution, and the scalar contribution, respectively. The upper right plot shows 
the same information and the cross-term (dot-dash line) with a limited y-axis to 
allow more detail to be seen. 



Fig. 4. The mKTr spectrum in data (error bars) comparing to the NR model with 
free Rr*(892)^ parameters (solid histogram) and the NR model with K*{892)^ pa¬ 
rameters fixed to the PDG [19] values (dashed histogram). 
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Table 2 

Summary of results on i^*(892)'^ parameters and contributions from non-(892)'^ 
sources in the decay —> K~Tr ~^obtained from the NR model. Fit result 

is compared to the current world averages and to the model with only Rr*(892)^. 
Limits on iLQ(1430)^ and iL*(1680)'^ contributions account for unseen decay modes. 



A* (892)° only 

FOCUS result 

Current values 

(892)0 (MeV/c^) 

895.61 ± 0.32 

895.41 ± 0.32;°;°! 

896.10 ± 0.27 [19] 

hic* ( 892)0 (MeV/c^) 

50.26 ± 0.81 

47.79 ± 0.86 + 

50.70 ± 0.60 [19] 

ro (GeV-i) 

14.1 ± 5.7 

3.96 ± 0.54 + |^ji^ 

3.40 ± 0.67 [17] 

r(D+^/s'-7r+M+DNR 

T(D+ ->K-■K+^l+v) ^ 

T{D+ 

V{D+ ^K--K+^l+v) 
r(D+^¥(5 (1430)0^+i^) 
V{D+ ^K--K+^l+v) 


5.30 ± 0.74 + °;°° 

< 4.0% @ 90% CL 

< 0.64% @ 90% CL 

~5[7] 

8.3 ± 2.9 [4] 

Confidence level (%) 

0.21 

66.0 



average value. We obtain a Blatt-Weisskopf parameter consistent with the 
value obtained by LASS [17]. We also limit possible additional Kn resonances 
present in semileptonic decays. Our angular distribution is 

consistent with the effective-range scalar non-resonant phase shift obtained by 
LASS [17] as expected by Watson’s Theorem given the absence of other hnal 
state interactions. 
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